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ABSTRACT 


An  expression  for  the  ray  curvature  of  a wave  packet  is  derived  which  is 
suitable  for  use  in  wave  prediction  programs.  The  ray  curvature  of  the  wave 
packet  is  found  to  vanish  if  the  packet  direction  becomes  either  perpendicular 
or  parallel  to  the  wave  speed  contours,  assuming  the  wavelet  direction  is  not 
parallel  to  the  contours.  For  parallel  water  depth  contours,  this  means  that  as 
a bydron  moves  into  shoaling  water  refraction  tends  eventually  to  turn  the  hydron 
so  that  it  is  directed  either  perpendicular  of  parallel  to  the  shoreline.  The 
first  case  is  similar  to  monochromatic  waves.  It  is  the  result  for  hydrons  which 
begin  in  deep  water  if  the  angle  of  incidence  is  between  0°  and  74.8°  with  respect 
to  the  contour  normal.  However,  for  deep  water  angles  of  incidence  equal  to  or 
greater  than  74.8°  the  hydrons  are  turned  and  move  parallel  to  shore  in  water  of 
intermediate  depth.  The  packet  ray  curvature  approaches  infinity  as  the  wavelet 
direction,  but  not  the  hydron  direction,  becomes  parallel  to  the  wave  speed  contours. 
The  result  is  total  reflection  of  the  waves.  Total  reflection  occurs  if  a hydron 
is  moving  into  deeper  water  and  its  initial  direction  exceeds  a critical  angle. 

At  the  reflection  point  the  hydron  direction  becomes  perpendicular  to  the  water 
depth  contours.  Further,  the  hydron  velocity  goes  to  zero,  which  is  consistent  with 


a particle  concept, 
countered. 
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Introduction 


The  conventional  definition  of  group  speed  U has  been  defined  as 

U - dw/dk  (1) 

where  u>  in  the  angular  frequency  and  k Is  the  wave  number.  This  equation  de- 

fines the  speed  of  the  group  in  the  direction  of  the  wavelet  velocity.  The 
geometric  group  speed  G was  defined  by  Breeding  (1978)  as 

G - U cos  $ (2) 

where 

♦ ■ 0 - Y (3) 

The  direction  of  movement  of  the  wave  packet  is  represented  by  8,  and  the 
direction  the  wavelets  move  within  the  packet  is  specified  by  y. 

Breeding  (1978)  has  shown  that  wave  packets  refract  according  to  Snell's 
law  with  the  geometric  group  velocity.  This  refraction  law  determines  the  wave 
path  of  constructive  interference.  At  each  point  of  the  wave  packet  trajectory 
the  wavelet  direction  is  determined  by  Snell's  law  with  phase  velocity. 

In  this  paper  an  expression  is  derived  for  the  ray  curvature  of  a wave 
packet  and  its  properties  are  described.  Examples  of  gravity  water  waves  will 
be  presented  to  show  how  G and  U differ  and  to  demonstrate  the  important  features 
of  wave  packet  refraction.  In  addition,  wave  packet  and  monochromatic  trajectories 
will  be  compared  for  waves  which  begin  in  deep  water. 

Ray  Curvature  for  Wave  Packets 

The  ray  curvature  of  a ray  moving  with  phase  speed  v was  derived  by  Munk 
and  Arthur  (1952)  and  Arthur,  et  al  (1952)  as 


dY  1 , . 3v  3v. 

si  ■ » (8ln  T a - c°* T 17) 


where  x,y  are  the  Cartesian  coordinates,  y is  the  direction  of 


respect  to  the  positive  x-axis,  and  sy  is  the  arc  length  along 


(4) 

the  ray  with 
the  ray. 
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The  ray  curvature  <G  for  the  trajectory  of  a wave  packet  is  given  by 

d0  1 , . n 3G  . 9G  . ,,v 

k„  = — * — (sin  G r cos  G r—  ) (5) 

G dsG  G 3x  3y 

where  0 is  the  direction  of  the  packet  ray  with  respect  to  the  positive  x-axis 
and  sG  is  the  arc  length  along  the  ray. 

Packet  Ray  Curvature  for  Locally  Parallel  Wave  Speed  Contours 

In  deriving  an  expression  for  the  packet  ray  curvature  the  wave  speed  con- 
tours are  assumed  to  be  parallel.  This  assumption  greatly  simplifies  the  derivation. 
Further,  it  leads  to  an  expression  suitable  for  making  calculations  from  which  the 
important  properties  of  the  packet  ray  curvature  can  be  determined.  Accordingly, 
the  xy-coordinate  system  is  chosen  such  that  the  y-derivatives  of  v,  U,  and  G are 
zero.  The  space  derivative  of  G is  then 

dG  dl]  . „ . . , d0  dv  x ,,, 

Expressions  for  the  space  derivatives  of  0 and  Y are  obtained  from  (4)  and  (5) 
where 

dx  * ds  cos  y (7) 

v 

dx  » ds„  cos  0 (8) 

Thus  dY/dx  = [(dv/dx)/v]  tan  y (9) 

d0/dx  ■ [(dG/dx)/G]  tan  0 (10) 

After  (9),  (10),  and  (6)  are  substituted  into  (5)  and  the  result  is  simplified, 
the  packet  ray  curvature  is  found  to  be 

[(dU/dx)/u]  4-  [(dv/dx)/v]  tan  tan  y 
KG  esc  0 + tan  $ sec  0 (11) 

This  expression  for  the  packet  ray  curvature  can  be  used  in  wave  prediction 

programs  even  if  the  wave  speed  contours  are  not  everywhere  parallel.  From  a 

practical  viewpoint,  it  is  only  necessary  that  the  wave  speed  contours  be  locally 

parallel  in  the  vicinity  of  each  ray  point.  Eq.  (11)  is  used  by  Breeding  et  al  1978 


L 
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in  computing  wave  packet  trajectories  of  gravity  water  waves.  To  use  (11)  the 
coordinate  system  is  rotated  at  each  ray  point  so  that  the  y-axis  is  parallel 
with  the  water  depth  contours  and  the  positive  x-axis  is  directed  toward  deeper 
water. 

Properties  of  the  Packet  Ray  Curvature  and  Snell's  Law 

The  ray  curvature  of  a wave  packet  defined  by  (11)  exhibits  some  very  re- 
markable properties.  It  is  assumed  that  v,  U,  and  their  derivatives  are  contin- 

uous and  finite.  However,  under  various  conditions  the  trigonometric  terms  of 
the  equation  can  become  infinite  or  have  indeterminate  forms.  The  value  of  <G 
approaches  zero  as  the  wave  packet  direction  6 becomes  either  parallel  or  perpen- 
dicular to  the  wave  speed  contours,  provided  the  wavelet  direction  y is  not 
parallel  to  the  contours.  This  means  that  given  a sufficiently  long  path,  re- 
fraction tends  to  turn  the  wave  packet  so  that  it  is  directed  either  parallel  or 

perpendicular  to  the  wave  speed  contours.  If  6 is  neither  parallel  nor  perpendicular 
to  the  wave  speed  contours,  then  *G  approaches  infinity  as  y becomes  parallel  to 
the  wave  speed  contours.  In  this  case,  due  to  the  value  of  y,  the  wave  packet 
undergoes  total  reflection. 

To  determine  the  value  of  k when  there  are  indeterminate  forms  it  is  necessary 

u 

to  consider  the  variations  of  0 and  y as  the  indeterminate  forms  are  approached. 

For  example,  (11)  contains  an  indeterminate  form  when  0 becomes  perpendicular  to 
the  wave  speed  contours  while  y becomes  parallel  to  the  contours.  If  y approaches 
parallelism  to  the  contours  faster  than  0 approaches  the  perpendicular  to  the 
contours  the  value  of  ic  becomes  infinite. 

U 

The  relationship  between  0 and  y due  to  refraction  is  clearly  seen  by  inte- 
grating the  ray  curvature  expressions  (4)  and  (5).  The  y-derivatives  being  zero, 
integration  of  (5)  leads  to 

(sin  0)/[u  cos  (0-y)]  ■ C (12) 

which  is  Snell's  law  for  a wave  packet  where  C is  a constant.  Snell's  law  with 
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phase  velocity,  which  determines  y.  is  obtained  by  integrating  (4).  The  cosine 
term  in  (12)  can  be  replaced  by  the  identity  for  the  difference  of  two  angles 
and  the  terms  rearranged  to  yield 

tan  0 = (UC  cos  y)/(l  - UC  sin  y)  (13) 

where  the  variation  of  0 appears  only  on  the  left  side  of  the  equation. 

It  is  interesting  to  note  that  0 becomes  zero  if  y * (2m  + 1)  (u/2)  where 
m is  an  integer.  Thus  if  the  wavelet  direction  becomes  parallel  to  the  wave 
speed  contours  the  wave  packet  direction  becomes  perpendicular  to  the  contours . 
Further,  note  that  0 = (2m  + 1)  ( tt/2)  if  UC  sin  y = 1.  For  this  case  the  wave 
packet  direction  is  parallel  to  the  wave  speed  contours. 

Snell's  law  can  be  used  to  derive  an  expression  for  cos  (0-y) . Eq.  (13)  is 
substituted  into  the  identity  for  tan  (0-y)  and  the  result  is  simplified  to  obtain 

tan  (0-y)  = (UC  - sin  y)/cos  y (14) 

In  terms  of  initial  values,  Snell's  law  with  phase  velocity  can  be  written 

sin  y * v^  sin  y^  (15) 

where  v^  = (v/v^)  and  the  subscript  i denotes  an  initial  value.  Before  refraction 
it  is  assumed  that  0^  = y^.  Then 


C = (Sin  Yi)/Ui 


(16) 


When  (15)  and  (16)  are  substituted  into  (14),  it  is  found  that 

(Ur  - vr)  sin  Yi 


tan  (0-y) 


2 2^ 
(1  - vf  sin  y±) 


(17) 


where  Uf  = (U/U^) . This  result  can  be  transformed  by  the  use  of  an  identity  to 


2 2 

(U  - v ) sin  y. 

cos  (0-y)  *>  [l+  — ^2 5 “ 1 

1 - vr  sin  Yi 


(18) 


Eqs.  (1),  (2),  and  (18)  provide  a useful  means  of  computing  the  geometric  group 
speed. 
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Examples  of  Wave  Packet  Refraction  - Hydrons 

To  demonstrate  the  properties  of  wave  packet  refraction,  examples  of  gravity 
water  waves  will  be  considered.  Gravity  water  waves  are  particularly  suited  as 
examples  because  of  their  highly  dispersive  nature.  The  term  'hydron,'  suggested 
by  Purser  and  Synge  (1962)  and  Synge  (1962),  will  be  used  to  denote  the  wave 
packet  of  water  waves. 

Waves  from  Deep  Water 

In  Figure  1 hydron  trajectories  are  shown  for  waves  beginning  in  deep  water 
(water  depth  greater  than  one  half  the  wavelength) . The  water  depth  contours 
are  parallel.  Initially  6^  * where  each  initial  direction  indicated  on  the 
figure  is  the  angle  between  the  hydron  velocity  vector  and  the  normal  to  the 
depth  contours.  Regardless  of  the  wave  period,  for  deep  water  angles  of  incidence 
between  0°  and  74.8°  the  hydrons  follow  paths  such  that  the  angles  increase  to  the 
depth  of  the  geometric  group  speed  maximum  (see  Figure  3),  then  undergo  a point 
of  inflection,  and  then  decrease  shoreward.  As  a hydron  approaches  shore  its 
direction  becomes  perpendicular  to  the  wave  speed  contours  and  the  packet  ray 
curvature  approaches  zero.  For  deep  water  angles  of  incidence  equal  to  or  greater 
than  74.8°  the  hydron  trajectories  turn  and  move  parallel  to  shore  in  water  of 
intermediate  depth.  As  the  hydron  direction  becomes  parallel  to  shore  the  packet 
ray  curvature  tends  to  vanish;  this  is  apparent  in  ray  number  5. 

For  comparison,  monochromatic  rays  are  shown  in  Figure  2 for  the  same  conditions 
considered  in  Figure  1.  For  large  incident  angles  there  is  a striking  difference 
between  hydron  and  monochromatic  trajectories.  Whereas  all  the  hydron  rays  do 
not  reach  shore  all  the  monochromatic  rays  do.  Note  that  the  wavelet  direction 
at  each  point  along  a hydron  path  in  Figure  1 is  the  same  as  the  direction  of  the 
corresponding  monochromatic  ray  at  the  corresponding  water  depth. 

It  Is  interesting  to  compare  the  values  of  G and  U when  they  differ  due  to 
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Figure  1.  Hydron  trajectories  for  a 20  sec  wave  period  for  waves  beginning 
in  deep  water.  The  water  depth  contours  are  parallel,  the  scale  of  the  plot 
is  1 cm  = 4.87  km,  and  the  sounding  water  depths  are  in  meters.  The  initial 
hydron  direction  is  shown  for  each  ray  and  is  the  angle  between  the  hydron 
velocity  vector  and  the  normal  to  the  depth  contours. 
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refraction.  As  an  example,  for  gravity  water  waves  which  begin  In  deep  water  it 
is  found  that 


U 


tanh  kh 
2kh 


(1  + 


sinh  2kh  ^ vr 


(19) 

(20) 


where  h is  the  water  depth.  When  (19)  and  (20)  are  substituted  into  (18)  and 
the  result  is  substituted  into  (2)  it  is  found  that 


G - u[l  + 


(kh  sin  y^  sech^  kh)^ 

1 - (sin  y^  tanh  kh)^ 


1 


(21) 


The  ratio  of  the  geometric  group  speed  to  its  initial  deep  water  value  is 
presented  for  several  incident  angles  in  Figure  3.  The  initial  hydron  directions 
are  defined  as  in  Figure  1.  The  curve  for  y^  ■ 0°  is  the  same  as  obtained  for 
the  ratio  of  U to  the  value  of  U in  deep  water.  The  amount  by  which  the  other 
curves  differ  from  it  is  a measure  of  the  difference  between  G and  U. 

Inspection  of  Figure  3 shows  for  a given  y^  that  the  maximum  of  G/G^  occurs 
at  a greater  value  of  kh  than  does  the  minimum  of  G/U.  An  increase  in  y^  causes 
a shift  in  both  the  minimum  of  G/U  and  the  maximum  of  G/G^  to  larger  values  of 
kh.  Further,  the  maximum  peak  tends  to  get  flattened  out.  The  curve  for  y^ 

3 74.8°  is  seen  to  stop  abruptly  at  the  maximum  value  of  G/G^. 

When  y^  ■ 30°  the  maximum  percentage  difference  of  G from  U is  2.70%.  When 
y^  » 45°  the  value  is  5.91Z,  for  y^  ■ 60°  it  is  10.03Z,  and  for  y^  ■ 74.79°  the 
value  is  14.27Z. 

Reflection  Points 

To  obtain  a reflection  point  it  is  necessary  that  the  waves  propagate  into 
deeper  water  and  that  the  initial  direction  of  the  hydron  (0^  • y^)  exceed  a 
critical  angle.  The  reflection  point  occurs  at  an  intermediate  water  depth  when, 
through  refraction,  the  wavelets  are  turned  parallel  to  the  wave  speed  contours. 

In  Figure  4 two  rays  are  shown  in  which  the  wave  period  is  20  sec  and  the 
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Figure  3.  The  variation  of  the  geometric  group  speed  with  kh  for 
hydrons  beginning  in  deep  water.  The  initial  value  of  the  geometric 
group  speed  is  G . . The  hydron  directions  are  defined  as  in  Figure  1 


Figure  4.  Hydron  trajectories  for  a 20  sec  wave  period  for  waves  which 
begin  at  an  intermediate  water  depth.  The  scale  of  the  plot  is  1 cm  - 
3.10  km  and  the  sounding  water  depths  are  in  meters.  The  initial  hydron 
direction  is  shown  for  each  ray  and  is  defined  as  in  Figure  1. 
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initial  water  depth  is  15  m.  For  this  case  a reflection  point  occurs  if  0^  22.2°. 

Ray  number  1 reaches  deep  water  since  0^  ■ 15°.  For  ray  number  2,  0^  * 23°, 
and  a reflection  point  occurs  at  a water  depth  of  200  m.  The  variation  of  the 
hydron  and  wavelet  directions  with  water  depth  for  this  ray  are  shown  in  Figure 
5.  The  wavelet  angle  increases  continuously  with  water  depth  and  becomes  parallel 
to  the  water  depth  contours  at  the  reflection  point.  The  hydron  angle  first  in- 
creases, goes  through  a maximum,  and  then  becomes  perpendicular  to  the  water  depth 
contours  at  the  reflection  point.  This  is  clearly  seen  in  Figure  4. 

Plots  of  how  G and  U vary  with  water  depth  in  approaching  the  reflection 
point  are  shown  in  Figure  5.  Both  speeds  have  a maximum  value.  However,  there 
is  a marked  contrast  between  G and  U at  the  reflection  point  where  the  geometric 
group  speed  is  zero. 

In  Figure  6 the  ratio  of  the  packet  ray  curvature  to  its  initial  value  is 
sketched  as  a function  of  the  wavelet  direction  for  ray  number  2 in  Figure  4.  The 
most  striking  feature  of  this  curve  is  that  within  about  2°  of  the  reflection 
point  the  packet  ray  curvature  suddenly  goes  to  infinity. 

It  is  interesting  to  observe  that  the  velocity  of  the  wave  packet  goes  to 
zero  at  the  reflection  point.  This  is  exactly  what  the  velocity  of  a particle 
should  do  when  there  is  reflection.  Therefore,  in  water  waves,  as  in  quantum 
mechanics,  the  wave-particle  duality  is  encountered. 

Summary 

Refraction  causes  a wave  packet  (hydron)  trajectory  to  become  directed  either 
parallel  or  perpendicular  to  the  water  depth  contours.  In  either  case  the  packet 
ray  curvature  will  vanish.  For  hydrons  propagating  toward  deep  water,  if  the 
initial  direction  exceeds  a critical  angle  total  reflection  occurs.  At  the  re- 
flection point  the  wavelet  direction  becomes  parallel  to  the  wave  speed  contours, 
the  hydron  direction  becomes  perpendicular  to  the  contours,  the  geometric  group 
velocity  goes  to  zero,  and  the  packet  ray  curvature  becomes  infinite. 
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Figure  5.  Variations  of  the  wavelet  direction  y,  hydron  direction  9, 
the  speed  U "dw/dk,  and  the  geometric  group  speed  G ■ U cos  (0-y) 
up  to  the  reflection  point  for  ray  number  2 in  Figure  4. 
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SPEEDS  G AND  U (M/SEC) 


Figure  6.  The  ratio  of  the  hydron  ray  curvature  to  its  absolute 
initial  value  as  a function  of  the  wavelet  direction  up  to  the 
reflection  point  for  ray  number  2 in  Figure  4. 
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